Strain effects on ferroelectric polarization and magnetism in orthorhombic HoMnOs 
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Aiming at increasing the ferroelectric polarization in AFM-E ortho-HoMnOa, we investigate the 
in-plane strain effects on both the magnetic configuration and the polarization by means of density 
functional theory calculations and model Hamiltonian approaches. Our results show that the net 
polarization is largely enhanced under compressive strain, due to an increase of the electronic con- 
tribution to the polarization, whereas the ionic contribution is found to decrease. We identify the 
electron-lattice coupling, due to Jahn- Teller (JT) distortions, and its response to strain to be re- 
sponsible for the observed behavior. The JT-induced orbital ordering of occupied Mn-Cg electrons in 



alternating 3x — r /3y — r orbital states at equilibrium changes to a mixture with x 



states under in-plane compressive strain. The asymmetric hopping of Cg electrons between Mn ions 
along zig-zag spin chains (typical of the AFM-E spin configuration) is therefore enhanced under 
strain, explaining the large value of the polarization. We reproduce the change in the orbital order- 
ing pattern in a degenerate double-exchange model supplemented with electron-phonon interaction. 
In this picture, the orbital ordering change is related to a change of the Berry phase of the Cg 
electrons, which in turn causes an increase of the polarization, whose origin is purely electronic. 

PACS numbers: 75.85.-|-t, 75.50.Ee 
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I. INTRODUCTION 

Orthorhombic rare-earth manganites iiMnOa, which 
crystalize in the Pbnm structure, represent a proto- 
typical class of "improper ferroelectrics" ,— where elec- 
tric dipoles are induced by a frustrated magnetic or- 
dering. Among this family of compounds, TbMnOa 
and DyMnOa show non-collinear magnetic structures, 
where spin-orbit coupling (SOC) via the Dzyaloshinskii- 
Moriya (DM) interaction is responsible for the weak 
ferroelectric polarization (P < 0.1 /xC/cm^) and the 
magnetoelectricityJi^ On the other hand, the on- 
set of a ferroelectric state induced by an exchange- 
striction mechanism (such effect is often called "inverse 
Goodenough-Kanamori interaction" §■) was theoretically 
proposed for HoMnOa in the antiferromagnetic E-type 
(AFM-E) spin confIguration /'^'^°i^^ and later experimen- 
tally confirmed. In the framework of the degenerate 
double exchange model;^ it has been shown that a dis- 
placement mechanism, based on the difference - in terms 
of bond-angles and bond-lengths - between Mn atoms 
with parallel and antiparallel spins, takes place for in- 
equivalent bridge oxygens, leading to a net ferroelectric 
polarization (P) directed along the a axis. The mecha- 
nism was confirmed by a density functional theory (DFT) 
study,* which also reported a comparable contribution 
arising from a purely quantum mechanical mechanism, 
resulting in a large value of the total ferroelectric polar- 
ization, F sa 6 /iC/cm^. Further calculations based on 
the Heyd-Scuseria-Ernzerhof hybrid functional report a 
value of the ferroelectric polarization of « 2 /uC/cm^j^ 
Microscopically, the electronic contribution is explained 



in terms of the asymmetric hopping of Mn-e^ electrons 
between neighboring Mn ions, where the "one-way" di- 
rection of hopping is determined both by the orbital or- 
dering and by the E-type spin configurationi^S As a sec- 
ond effect, the ions move to enhance the polarization in- 
duced by the asymmetric hopping; in fact, the Mn-O-Mn 
angle between the ferromagnetically coupled Mn spins 
increases after ferroelectric ionic relaxations, resulting in 
an increase of the Cg hopping integrals. 



For an easier detection of the polarization in exper- 
iments, an increase of the polarization in this system 
is desired. A possible path to achieve this is by in- 
creasing the Mn-Cg hopping amplitude by modulating 
the (ab) in-plane Mn-O-Mn bond angle cj) (average value 
of 0=144° for HoMnOgi^), as well as the Mn-0 bond 
lengths. However, the enhancement of (f> may lead, ac- 
cording to the Goodenough-Kanamori rule,^^ to the sta- 
bilization of other magnetic phases, i.e. the non-collinear 
configuration, with a smaller P, observed in TbMnOa 
((/)=145°) and the paraelectric AFM-A phase in LaMnOs 
((/>=155°). Another possible recipe to increase the polar- 
ization P is to shorten the bond lengths in the AFM-E 
magnetic ground state, hopefully keeping the angle un- 
changed in such a way that the AFM-E phase does not 
turn unstable. In this context, here we theoretically in- 
vestigate strain effects both on the magnetic configura- 
tion and on the polarization for orthorhombic HoMnOs, 
where the in-plane lattice constants a and b are artifi- 
cially modified in order to mimic the compressive/tensile 
modes acting in thin films grown on different substrates. 
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II. STRUCTURAL AND COMPUTATIONAL 
DETAILS 

Orthorhombic HoMnOscrystalizes in the Pbnm struc- 
ture (corresponding to Pnma in the standard setting; in 
this paper, we chose c as the longest axis). Here, strong 
distortions with respect to the ideal cubic perovskite, 
i.e. GdFeOa type tilting of MnOg octahedral cages and 
Jahn- Teller (JT) distortions, stabilize the orbital order- 
ing of Mn-Bg states, where the orbital order is expressed 
as (3a;^ — r^)/(3y^ — r^). Although the hexagonal non- 
perovskite phase is more stable than the orthorhombic 
phase in HoMnOa, the transition to the metastable or- 
thorhombic phase can be obtained by high-pressure syn- 
thesis. 

DFT simulations were performed using the VASP 
codeii and the PAW pseudopotential&l^^ within the GGA 
formalism. (Note that for ortho-i^MnOa, it is already 
known that the GGA+U treatment^^ worsens the agree- 
ment with experimental results compared to bare GGA, 
as regards the optimized atomic structure and the mag- 
netic stabilizationiiSi^) The cut-off energy for the plane- 
wave expansion of the wave-functions was set to 400 eV 
and a k-point shell of (4, 2, 3) was used for the Brillouin 
zone integration according to the Monkhorst-Pack special 
point scheme. The in-plane lattice constants correspond- 
ing to the unstrained structure are taken as a = 5.2572 A 
and b = 5.8354 A, equal to the experimental values.l-. In 
order to mimic strain, these values were progressively 
changed from -4% (compressive strain) to 4% (tensile 
strain). The lattice constant c was optimized for each 
(a,b) value, as well as the internal coordinates. In the 
unstrained case, the theoretically optimized c is found to 
be equal to 7.3952 A, slightly underestimated compared 
to the experimental value of 7.3606 A. In the AFM-E 
magnetic configuration ( AFM q vector along b) , the unit 
cell is doubled along b (40 atoms/cell). Therefore, the 
crystal structure is optimized in the non-centrosymmetric 
P2imn symmetry, allowing for a ferroelectric polariza- 
tion along X (Px). (Without spin-orbit coupling, the 
magnetic group is not considered.) A schematic repre- 
sentation of the structure in the AFM-E magnetic con- 
figuration is shown in Fig. [T] The electronic polarization 
^'Borry=^tot is Calculated using the Berry phase method^ 
and compared to the classical point charge model (PCM) 
result, where PpcM=-Pion is calculated as the product 
of the nominal ionic valence and of the ionic displace- 
ments. In order to discuss the stability of different mag- 
netic configurations, the internal atomic coordinates were 
fully optimized in the AFM-A configuration, character- 
ized by in-plane ferromagnetic and inter-plane antiferro- 
magnetic couplings. Using this relaxed centrosymmetric 
structure, the total energies of several magnetic config- 
urations (antiferromagnetic AFM-A, AFM-C, AFM-G, 
AFM-E, AFM-E* and ferromagnetic FM configurations) 
were compared, and the exchange coupling integrals Jij 
were calculated by mapping the system onto a Heisenberg 
model (see Ref. J^for details). 




FIG. 1: (Color online) Spin density of the Mn-Cg states 
(within an energy range up to 1 eV below the Fermi energy) 
in AFM-E configuration without strain, projected in the (ab) 
plane, referred in the text as in-plane. Blue (dark) and red 
(light) isosurface show up- and down-spin state (at ±0.1 e~), 
respectively. Due to the exchange striction, the Mn-O-Mn 
bond angles 0^ between parallelly-coupled Mn-spins and 
between anti-parallelly-coupled Mn-spins are different, as well 
as the Mn-O bond distances F, and s^, s""^ . The direc- 
tion of the net polarization (P || —a) is denoted by a big red 
arrow. 



III. STRAIN EFFECT ON ATOMIC 
STRUCTURE 

Fig. [5] shows the dependence of the total energy on 
strain. The deviation from the origin of the energy mini- 
mum, giving an estimate of the error of bare GGA treat- 
ment with respect to experimental results, is less than 
1% (rather typical for DFT calculations). The band gap, 
given by the JT-splitting of Mn-e^ states, is also found 
to decrease under compressive strain in the (ab) plane. 

In Fig. 13] and Fig. 0] we show how the change in the 
a and b lattice constants affects the optimized value of 
the lattice constant c and the lattice volume, as well as 
the Mn-0 bond lengths and Mn-O-Mn bond angles in 
the AFM-E configuration. Under in-plane compressive 
strain, the c lattice constant is expanded, while the vol- 




ab strain (%) 

FIG. 2: (Color online) Calculated variation of the total energy 
difference and the band gap with optimized atomic structure 
in AFM-E configuration under in-plane strain. 
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FIG. 3: (Color online) Calculated variation of equilibrium c 
and lattice volume V under in-plane strain, a and b being 
fixed for a given strain. 



The inter-plane Mn-O-Mn bond angles, 0*°^ and 
show a rapid increase under compressive strain. Recall- 
ing the fact that the hopping integral between Mn-e^ 
states depends both on the bond angle (j> and on the 
bond length d — the overlap integral b is described as — 
b = cos((180° — 0)/2)/d'^-^ — , apparently oxygen ions are 
displaced in such a way that the hopping energy is max- 
imized. On the other hand, the difference [{/jP-c^^^Pj of the 
bond-angle between Mn ions with parallel and antipar- 
allel spins (proven to be relevant for the ionic contribu- 
tion to the polarization) slightly increases under tensile 
strain, probably because of the smaller elastic energy that 
binds oxygen atoms in a relatively larger space in the lat- 
tice. Such a reduction of the elastic energy is expected to 
enhance ionic displacement and, consequently, the ionic 
contribution to polarization, Pion- 
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FIG. 4: (Color online) Variation of (a) the Mn-O bond 
lengths, in-plane I, s (long axis and short axis in JT-distorted 
MnOs octahedron) and inter-plane m and (b) the Mn-O-Mn 
in-plane bond angles (j) and inter-plane bond angles (p in the 
AFM-E configuration under (ab) in-plane strain. 



ume is contracted. The Mn-0 bond lengths and the Mn- 
O-Mn bond angles depend on the relative alignment of 
the Mn spins, parallel or antiparallel, as a consequence 
of exchange striction. Remarkably, the in-plane bond an- 
gles, 0^ and 4)°'P, are kept rather constant under strain, 
whereas the bond lengths are much more modulated. 



IV. MAGNETIC CONFIGURATION 

In this section we focus on the magnetic stabilization 
under strain. As shown in Fig. [5j the AFM-E phase is 
stable under tensile strain, but, as a general trend, un- 
der compressive strain (at -4%) none of the AFM phases 
is any longer stable with respect to the ferromagnetic 
phase. This trend is consistent with the evolution of the 
band gap, which closes towards a metallic state under 
compressive strain. Specifically, the double exchange in- 
teraction between the Mn-e^ states is enhanced in this re- 
gion of strain. In fact, the in-plane nearest-neighbor FM 
coupling J|"" progressively increases under compressive 
strain. The AFM-E configuration, gaining energy both 
from the double exchange interaction along the zigzag 
spin chains and from the ferroelectric ionic displace- 
ments, is more stable than the FM and AFM-A phases 
within a region of strain ranging from -3% to 4%. It is 
also worth noticing that the inter-plane next-neighbor 
coupling J]_ rapidly changes sign under compressive 
strain, turning from AFM to FM. The reason for the 
emergence of a FM inter-plane coupling will be explained 
in the next section. From the calculated exchange inter- 
actions, Jij, we can also estimate the critical temperature 
Tc in the mean-field approximation for the Heisenberg 
model: we find that the Neel temperature for the AFM- 
E phase is ks Tn = 2/3{jp + J^) and that the Curie 
temperature is ks Tc = 2/3{2Jf + + Jl + 4 ) for 
the FM phase, with ks the Boltzmann constant. Their 
dependence on strain is shown in Fig. [S] 



V. FERROELECTRICITY AND JAHN-TELLER 
DISTORTION 

Finally, we consider the strain effect on the ferroelec- 
tric polarization in the AFM-E phase. In order to discuss 
the evolution of P as a function of strain, we calculate 
the polarization also at -4% and -3% strain, where the 
AFM-E phase is not magnetically stable with respect to 
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FIG. 5: (Color online) Variation of (a) magnetic stabiliza- 
tion (total energy of AFM-A, C, G, E and E* phases with 
respect to FM in the atomic structure fixed to the opti- 
mized AFM-A phase, and the total energy of the AFM-E 
phase for which the atomic structure was optimized in the 
non-centrosymmetric phase), (b) Interatomic exchange in- 
teraction (Jij); The first nearest-neighbor in-plane J^^'^, next 
nearest-neighbor in- plane ^i""", first nearest- neighbor inter- 
plane Jx, next nearest-neighbor inter-plane J^. Positive 
(negative) sign denotes AFM (FM) coupling, (c) Ordering 
temperatures of the FM and AFM-E phases within the mean- 
field approximation. 




ab strain (%) 

FIG. 6: (Golor online) Variation of the ferroelectric polariza- 
tion, Pscrry (elcctronic+iouic contribution) and PpcM (ionic 
contribution), with strain. 



the FM one. As shown in Fig. [51 interestingly, it turns 
out that -P(=PBorry) IS enhanced at both limits of strain: 
i.e. under tensile strain, the contribution to P arising 
from ionic displacements (=Ppcm) increases, whereas the 
electronic contribution to P {=PBerTy-PpCM) decreases, 
with an opposite trend in the opposite limit of compres- 
sive strain. The former trend is easily understood: in a 
larger space, oxygen ions are more freely displaced by the 
exchange striction, resulting in an enhancement of Pion, 
as discussed in Sec. IIIIl The latter is less easy to under- 
stand intuitively. One possibility is that the hopping of 
Mn-Cg electron increases due to the shorter bond length, 
so that the asymmetric hopping is increased, thus ex- 
plaining the enhanced electronic contribution to the net 
polarization. In Fig. [71 the bandwidth w of Mn-e^ states, 
which reflects the hopping integral between Mn-Cg states 
at neighboring Mn-sites, is shown to be unchanged re- 
gardless of strain. However, the JT splitting is changed 
significantly, as one may have expected by looking at the 
trend of Mn-0 bond lengths shown in Fig. [H where the 
long axis I and short axis s get closer under compres- 
sive strain. Besides, at —4% of strain, I and m becomes 
similar, and as a consequence changes sign from neg- 
ative to positive, as shown in Fig. [SI This change of 
the JT distortion is reflected in the shape of the ordered 
orbitals. Figure [SI shows the change of the orbital shape 
under compressive strain. The — orbital acquires 
additional character of the x^—z^ orbital (in a local frame 
X = I, y = s, z = m). As a consequence of the orbital 
change, an inter-plane Mn-0 bonding through Cg elec- 
trons is formed (c.f. small weight of the charge density 
on the apical oxygen atom in Fig. [H^a) but not in (b)); 
this might also be related to the fact that inter-plane 
magnetic coupling turns FM under compressive strain. 
Regarding the enhancement of the electronic contribu- 
tion to P, we should consider the charge transfer caused 
by Mn electrons at occupied eg state, which hop into 
unoccupied eg state at neighboring Mn site. It is ex- 
pected that the charge transfer is enhanced when the 
JT-splitting is small. 
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FIG. 7: Partial density of states (DOS) at Mn site under no 
strain and -4% strain. JT, w and CEF denote JT splitting at 
Mn-eg state, band width (oc hopping integral t) of Cg occupied 
state, and crystal field splitting between Cg and t2g states. 

VI. MODEL HAMILTONIAN 

The change in orbital ordering as a response to com- 
pressive strain found in the DFT framework can be ac- 
counted for in a model Hamiltonian approach. We con- 
sider a degenerate double exchange model, where infi- 
nite Hund coupling between the t2g and the Cg electrons 
is assumed.^'' This allows us to neglect spin degrees of 
freedom of conduction electrons, which arc fixed by the 
underlying magnetic structure of the core spins. In the 
E-type configuration, electrons can move along zig-zag 
chains on a background of parallel (ferromagnetically 
coupled) core spins, which experience an antiferromag- 
netic interchain coupling. Orbital ordering can then be 
stabilized by the Jahn- Teller coupling to the lattice. The 
Hamiltonian reads as: 

^ = - E cos e., cl c^. ^ -f J ^ s, . s, 

+ i?jT^ [2(92 + qnTzz) + qli + 93 J • (1) 

i 

We adopt the standard notation, with (c^^) creation 
(annihilation) operators for electrons in the bands a = 




ab strain (%) 

FIG. 8: (Color online) Variation of the Jahn- Teller vector Q 
= (Q2, Q3)= [-^{l ~ s), ^(2m - / - s)y^ in the AFM-E 
configuration, where I, m and s are averaged, under in-plane 
strain. 




(a) ab -4% (b) ab +2% 



FIG. 9: (Color online) Isosurface of charge density of occupied 
Mn-Cg state (within an energy range up to leV below Fermi 
energy) under (a) -4% and (b) -1-2% strain. 

a, b (stemming respectively from d^^-y^ , d^z^~r^ orbitals 
on Mn'^"'"); t"^^ are the oxygen-mediated hopping ampli- 
tudes between nearest-neighbors manganese atoms, i.e. 
t'ia = tla = 3to/4, ft, = tl = io/4, ft, = ~tl^^ = V3to/4 
with to — (pda)^ as the energy reference.^^ The &ij are 
the angles between neighboring Mn-core spins 5*^ and Sj 
and J is the exchange coupling which they experience. 
The JT modes are expressed through the dimensionless 
121, Qsi (= (k/g) Qi) and Ejt — g^f^k is the static Jahn- 
Teller energy; = Y^ap 4a '^a/s '^i,i3 are the orbital pseu- 
dospins, where ct^^ are the Pauli matrices. 

We evaluated the optimal q2i , q^i keeping the t2g spins 
in E-type configuration, and we found that in the un- 
doped case (one electron per Mn site) an orbital order- 
ing occurs as soon as the JT coupling is large enough 
to induce a finite distortion. The character of orbitals 
involved in the ordered phase can be visualized by eval- 
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FIG. 10: (Color online) Jahn- Teller phase at neighbor- 
ing sites as a function of strain as obtained in DFT (left), 
compared with orbital phase 6i evaluated in model Hamilto- 
nian approach as a function of crystal field with reversed sign 
(right). Dotted lines at 7r/3 (2n/3) and 5n/3 (47r/3) corre- 
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and y 



uating the average occupancy of the atomic orbitals, de- 
fined as pai = (0L<?^ai)' where (j>ai = - s'm{0i/2) Cai + 
cos{9i/2) Cbi and a labels each of the possible e^-electron 
orbitals with the corresponding angle 6i (e.g. Oi = 
27r/3(47r/3) corresponds to 3a:^ - (3y^ - r^) orbital, 
see for instance Ref. [2^ . 

The effect of the in-plane strain can be mimicked by 
tlie introduction of a crystal-field term -|-A r^j in the 
Hamiltonian ([1]); if we consider a compression in the ab 
plane, then A should be chosen positive. Another ef- 
fect of compression is of course to increase the elastic 
stiffness of MnOg octahedra, therefore one should expect 
a decrease of JT static energy (as indeed observed in 
the reduction of the JT gap found in DFT calculations). 
We found that, upon increasing A, the orbital ordering 
is not affected in its alternating pattern, but there is a 
change in the weight of involved atomic orbitals, namely 
the weight of 
3a;2 - (3?/^ - 
the Cg-electron bands is reduced even with fixed JT cou- 
pling. To visualize the evolution of orbitals involved in 
the ordered pattern as a function of strain we can evalu- 
ate the angle Ot = tan~^((Ta;i)/(Tzi)) and compare it with 
the Jahn- Teller phase — tan~^((52i/Q3i) extrapolated 
from DFT calculations. As shown in Fig. [TUl we find a 
good qualitative agreement between the two phases. 

We have also checked our results by considering the ef- 
fect of cooperative distortions of oxygens. This has been 
done by optimizing the Hamiltonian ([Ij with respect to 
oxygen displacements {u} from their equilibrium posi- 
tions rather than optimizing the local distortions of the 
MnOg octahedron labeled by {q}: in this case the strain 
can be mimicked by changing the equilibrium positions 
of oxygens, specifically decreasing the MnO distances in 



{y^ — z^) is increased with respect to 
At the same time, the gap between 



Oh 




FIG. 11: (Color online) Electronic polarization as a function 
of -A (A is the crystal field). 



the ah plane with respect to the equilibrium position of 
oxygens along the c axis. This approach provides results 
which are consistent with those obtained in the nonco- 
operative treatment of Jahn- Teller effect. 

We focus then on possible electronic polarization. Due 
to the AFM-E magnetic configuration, at T = elec- 
trons can only hop between sites with ferromagnetically 
aligned core-spins. This implies that hopping processes 
can occur only al ong one-dimensional zig-zag chains. As 
discussed in Ref. [2a, the presence of Jahn- Teller centers 
in this one-dimensional system can induce a Berry phase 
in the Bloch function, intimately related to the orbital or- 
dering induced by JT distortions, from which electronic 
polarization can be inferred. We follow the prescription 
described in Ref. 27 and evaluate the polarization within 

the chains as Pei = lim — lmln(4'o| e*~'^|^'o), where 
L— >oo 27r 

L is the length of the one-dimensional chain (thermody- 
namic limit is taken at the end of the calculation), is 
the ground-state wavefunction of model ([T]) with opti- 
mized q2i,q3i, and X — in is the electronic position 
operator. The onset of the orbital ordering within the 
zigzag chains induces a finite polarization whose origin 
is fully electronic, the ions position being fixed in the 
system.— As shown in Fig. [TTl we find that Pel is in- 
deed enhanced by increasing A (i.e. as a consequence of 
a compressive strain in our model approach) , in excellent 
qualitative agreement with the DFT results. 



VII. CONCLUSIONS 

The ferroelectric and magnetic properties of or- 
thorhombic HoMnOs were investigated by a combined 
density functional theory and model Hamiltonian study. 
We found an increase in the electric polarization when 
the system is subjected to in-plane compressive strain. 
This is closely related to the modifications in the elec- 
tronic structure, namely the orbital ordered states — 

/Sy^—r"^ gain additional character of the — / y'^ — 
states. Beside leading to an increase of the polarization, 
the orbital mixing infiuences the magnetic stability; thus, 
a crossover from the AFM-E to the FM phase occurs 
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around -3% compressive strain. 

Acknowledgments 

The research leading to these results has received fund- 
ing from the European Research Council under the EU 
Seventh Framework Programme (FP7/2007-2013) / ERG 
grant agreement n. 203523. Computational support from 



the Swedish National Infrastructure for Scientic Com- 
puting (SNIC), Caspur Supercomputing Center (Rome) 
and Cineca Supercomputing Center (Bologna) is grate- 
fully acknowledged. B. S. acknowledges support from 
VR. O. E. acknowledges support from VR, SSF, ERC, 
and the European Comission, project number FP7-NMP- 
2008-EU-India2-233513. The charge density plots were 
performed using the program VESTAi^ 



* |diana.iusan@spin.cnr.it| 

^ S. W. Cheong and M. Mostovoy, Nat. Mater. 6, 13 (2007). 
^ A. P. Levanyuk and D. G. Sannikov, Sov. Phys. Usp. 17, 
199 (1974). 

^ M. lizumi, J. D. Axe, G. Shirane and K. Shimaoka, Phys. 
Rev. B 15, 4392 (1977). 

* M. Kenzelmann, A. B. Harris, S. Jonas, C. Broholm, J. 
Schefer, S. B. Kim, C. L. Zhang, S.-W. Cheong, O. P. 
Vajk, and J. W.Lynn, Phys. Rev. Lett. 95, 087206 (2005). 

^ I. A. Sergienko and E. Dagotto, Phys. Rev. B 73, 094434 
(2006). 

® K. Yamauchi, S. Picozzi, J. Phys. : Condens. Matter 21, 
064203 (2009). 

L A. Sergienko, C. Sen, and E. Dagotto, Phys. Rev. Lett. 
97, 227204 (2006). 

* S. Picozzi, I. A. Sergienko, K. Yamauchi, B. Sanyal and E. 
Dagotto, Phys. Rev. Lett. 77, 227201 (2007). 

^ A. Stroppa and S. Picozzi, Phys. Chem. Chem. Phys. 12, 
5405 (2010). 

K. Yamauchi, F. Freimuth, S. Blugel, and S. Picozzi, Phys. 
Rev. B 78, 014403 (2008). 

S. Picozzi and C. Ederer, J. Phys.: Condens. Matter 21, 
303201 (2009). 

^2 B. Lorenz, Y. Q. Wang, and C. W. Chu, Phys. Rev. B 76, 
104405 (2007). 

B. Lorenz, Y. Q. Wang, Y. Y. Sun, and C. W. Chu, Phys. 
Rev. B 70, 212412 (2004). 

V. Yu. Pomjakushin, M. Kenzelmann, A. Donni, A. B. 



Harris, T. Nakajima, S. Mitsuda, M. Tachibana, L. Keller, 
J. Mesot, H. Kitazawa, and E. Takayama-Muromachi, New 
J. Phys. 11, 043019 (2009). 

J. A. Alonso, M. J. Martnez-Lope, M. T. Casais, and M. 

T. Fernandez-Daz, Inorg. Chem. 39, 917 (2000). 
" J. B. Goodenough, Phys. Rev. 100, 564 (1955); J. 

Kanamori, J. Phys. Chem. Solids 10, 87 (1959). 
" G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 

(1996). 

^® P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

V.I. Anisimov, F. Aryasetiawan and A.I. Lichtenstein, J. 
Phys.: Cond. Mat. 9, 767 (1997). 

S. Picozzi, K. Yamauchi, G. Bihlmayer, and S. Bliigel, 
Phys. Rev. B 74, 094402 (2006). 
2^ R.D. King-Smith and D. Vanderbih, Phys. Rev. B 47, 1651 
(1993); R. Resta, Rev. Mod. Phys 66, 899 (1994). 
K. Momma and F. Izumi, J. Appl. Crystallogr., 41, 653 
(2008). 

J.-S. Zhou and J. B. Goodenough, Phys. Rev. Lett. 96, 
247202 (2006). 

E. Dagotto, T. Hotta, and A. Moreo, Physics Reports 344, 
1 (2001). 

J. C. Slater and G. F. Koster, Phys. Rev. 94, 1498 (1954). 
'^^ H. Koizumi, T. Hotta and Y. Takada, Phys. Rev. Lett. 80, 
4518 (1998). 

R. Resta, Phys. Rev. Lett. 80, 1800 (1998). 

P. Barone, K. Yamauchi, and S. Picozzi, 

larXiv:lQ09.336 9vl. 



